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ABSTRACT 

The recent robust and homogeneous analysis of the world's supernova 
distance-redshift data, together with cosmic microwave background and baryon 
acoustic oscillation data, provides a powerful tool for constraining cosmological 
models. Here we examine particular classes of scalar field, modified gravity, and 
phenomenological models to assess whether they are consistent with observations 
even when their behavior deviates from the cosmological constant A. Some mod- 
els have tension with the data, while others survive only by approaching the 
cosmological constant, and a couple are statistically favored over ACDM. Dark 
energy described by two equation of state parameters has considerable phase 
space to avoid A and next generation data will be required to constrain such 
physics. 
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Introduction 



A decade after the discovery of the acceleration of the cosmic expansion (jPerlmutter et al. 



19991 : iRiess et aLlll998l ) we still understand little about the nature of the dark energy physics 
responsible. Improved data continues to show consistency with Einstein's cosmological con- 
stant A, and in terms of a constant equation of state, or pressure to density, r atio w, the best 
fit to the data is w = -0.969lo:o63(stat)lo:o66(sys), where A has w = -1 JKowalski et al. 



20081 ). However, the magnitude of A required and the coincidence for it becoming dominant 
so close to the present remain unexplained, and an abundance of motivated or unmotivated 
alternative models fills the literature. Using the latest, most robust data available we exam- 
ine the extent to which data really have settled on the cosmological constant. 

The vast array of models proposed for dark energy makes comparison of every model 
in the literature to the data a Sisyphean task. Here we select some dozen models with 
properties such as well defined physical variables, simplicity, or features of particular physical 
interest. These embody a diversity of physics, including scalar fields, phase transitions, 
modified gravity, symmetries, and geometric relations. While far from exhaustive, they 
provide roadmarks for how well we can say that current data have zoomed in on A as the 
solution. 

For such comparisons it is critical to employ robust data clearly interpretable within 
these "beyond A" cosmologies. Geometric probes from the Type la supernovae (SN) distance- 
redshift relation, cosmic microwave background (CMB) acoustic peak scale shift parameter, 
and baryon acoustic oscillations (BAO) angular scale serve this essential role. Equally im- 
portant is confidence in the error estimates, incorporating systematics as well as statistical 
uncertainties. This has been studied in detail in the recent unifi ed analysis of the wo rld's 
published heterogeneous SN data sets - the UnionOS compilation (IKowalski et al.ll2008l ). 



This SN compilation includes both the large data samples from the SNLS and ESSENCE 
survey, the compiled high redshift SNe observed with the Hubble Space Telescope, a new 
sample of nearby SNe, as well as several other, small data sets. All SNe have been analyzed in 
a uniform manner and have passed a number of quality criteria (such as having data available 
in two bands to measure a color, and sufficient lightcurve points to make a meaningful fit). 
The samples have been carefully tested for inconsistencies under a blinded protocol before 
combining them into a single final data set comprising 307 SNe, the basis for this analysis. 
In this work the SNe da ta will be combined w ith the constraints obtained from the baryon 



acoustic oscillation scale (lEisenstein et a. 



and ground based CMB measurements (IKomatsu et al.ll2008l ). 



and from the fi ve year data release of WMAP 



In Section [2] we describe the general method for cosmological parameter estimation and 
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present a summary table of the various models considered and the statistics of the fit. 
Sections [314T2] then briefly describe the dark energy models, their parameters, and show the 
likelihood contours. The concluding discussion occurs in Section [T51 

2. Constraining Models 

Achieving informative constraints on the nature of dark energy requires restricting the 
degrees of freedom of the theory and the resulting degeneracies in the cosmological model 
being tested. One degree of freedom entering the model is the present matter density Qm- 
For the case of the spatially flat cosmological constant A model (or some of the other models 
considered below), this is the sole cosmological parameter determining the distances entering 
the supernova (SN) magnitude-redshift, baryon acoustic oscillation scale (BAO), and cosmic 
microwave background (CMB) shift parameter relations. 

Generally, further degrees of freedom to describe the nature of the dark energy, i.e. 
its equation of state (EOS), or pressure to density, ratio, are needed. In a few cases the 
EOS is parameter free, as in the A case where w = —1, or is determined by the matter 
density, as in some subcases below (such as the fiat DGP braneworld gravity model of [2]). 
One way to categorize models is by the number of independent EOS parameters, or general 
parameters beyond the matter density (so fiat A models have zero such parameters, A models 
with curvature have one). In general, current data can deliver reasonable constraints on one 
parameter descriptions of dark energy. 

In addition to exploring the nature of dark energy through its EOS, one might also 
include another parameter for the dark energy density, i.e. allow the possibility of nonzero 
spatial curvature. In this case individual probes then generally do a poor job constraining the 
model with current data, although the combined data from SN+CMB+BAO can sometimes 
still have leverage. Since crosschecks and testing consistency between probes is important 
(as particularly illustrated below in the DGP case), we consider spatial curvature only in 
the otherwise zero parameter cases of A and DGP, and for the constant EOS dark energy 
model. 

In the following sections we investigate various one parameter EOS models, discussing 
their physical motivation or lack thereof, and features of interest, and the observational 
constraints that can be placed upon them. In the last sections we also investigate some two 
parameter models of interest, with constrained physical behaviors and particular motivations. 
As a preview and summary of results. Table [1] lists the models, number of parameters, and 
goodness of fit for the present data. 
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The SN, CMB, and BAO data are combined by multiplying the likelihoods. Especially 
when testing models deviating from the cosmological constant one must be careful to account 
for any shift of the CMB sound horizon arising from violation of high redshift matter domina- 
tion on the CMB and BAO scales; details are given in UnionOS. Note that some doubt exists 
on the use of the BAO constraints for cos mologies other than ACDM, or possibly constant 



w. (IDick et al.ll2006l: iRvdbeck et al.l 120071 ) since ACDM is assumed in several places in the 
Eisenstein et al.l ( l2005l ) analysis, e.g. computation of the correlation function from redshift 
space, nonlinear density corrections, structure formation and the matter power spectrum, 
and color and luminosity function evolution. Properly, a systematic uncertainty should be 
assigned to BAO to account for these effects; however, this requires a complex analysis from 
the original data and we show only the statistical error. At the current level of precision, 
simplified estimates show this does not strongly affect the results, but such systematics will 
need to be treated for future BAO data. All figures use the likelihood maximized over all 
relevant parameters besides those plotted, and contours are at the 68.3%, 95.4%, and 99.7% 
confidence level. 

It is particularly important to note the treatment of systematic errors, included only 
for SN. We employ the prescription of UnionOS for propagation of systematic errors. This 
introduces a new distance modulus fi^^^ = fi + AMj + AM, which is simply the usual 
distance modulus fi = 5\og{HodL{z)), where ^^(z) is the luminosity distance and Hq the 
Hubble constant, shifted by a sample dependent magnitude offset AMj and a single sample 
independent magnitude offset AM added only for the higher redshift SNe [z > 0.2). The 
magnitude offsets AMj reflect possible heterogeneity among the SNe samples while the AM 
step from SNe at z < 0.2 to z > 0.2 allows a possible common systematic error in the 
comparison of low vs. high redshift SNe. Treating AMj and AM as additional fit parameters, 
one defines Xsys = X^+X]j(^^i/'^Mj^ + (AM/crA/)^ to absorb the uncertainty in the nuisance 
parameters, and ctm, and obtain constraints on the desired physical fit parameters 
that include systematic errors. This procedure of incorporating systematic errors provides 
robust quantification of whether or not a model is in conflict with the data and is essential 
for accurate physical interpretation. See UnionOS for further, detailed discussion of robust 
treatment of systematics within the current world heterogeneous SN data. 



3. Constant Equation of State 

Models with constant equation of state w within 20%, say, of the cosmological constant 
value w = —1, but not equal to —1, do not have much physical motivation. To achieve 
a constant equation of state requires fine tuning of the kinetic and potential energies of 
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Model 


Motivation 


Parameters 


X (Stat) 


X (sys) 


AL>um (natj 


gravity, zeropoint 




O 1 O 1 

olo.l 


OC\C\ C\ 

309. y 








Ax' (stat) 


Ax' (sys) 


ACDM 


gravity, zeropoint 




-1.1 


-1.3 


Constant w (flat) 


simple extension 


Qm, W 


-0.3 


-1.2 


Constant w 


simple extension 




-1.1 


-1.6 


Braneworld 


consistent gravity 




15.0 


2.7 


Doomsday 


simple extension 


^doom 


-0.1 


-0.7 


Mirage 


CMB distance 


Qrn, Wo 


-0.2 


-0.1 


Vacuum Metamorphosis 


induced gravity 




0.0 


0.0 


Geometric DE Riow 


kinematics 


ro,ri {flm,wo) 


0.1 


-1.1 


Geometric DE Rhigh 


matter era deviation 




-1.9 


-2.2 


PNGB 


naturalness 


^m, Wo, f 


-0.1 


-0.7 


Algebraic Thawing 


generic evolution 


Qrn, Wo, P 


-1.6 


-2.3 


Early DE 


fine tuning problem 


^m, Wo, 


-0.3 


-1.2 


Growing i/-mass 


coincidence problem 


Q Q 


-0.6 


-1.6 



Table 1: "Beyond A" dark energy models considered in this paper, together with ACDM 
models. Models are listed in the order of discussion, and the cosmological fitting param- 
eters shown. The x' of the matter plus cosmological constant case is given, and all other 
models fist the Ax' from that model. The values refer to the best fit to the joint data of 
SN+CMB+BAO; in the last column the SN systematics as analyzed in UnionOS are included. 
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a scalar field throughout its evolution. It is not clear that a constant u; 7^ —1 is a good 
approximation to any reasonable dynamical scalar field, where w varies, and certainly does 
not capture the key physics. However, since current data cannot discern EOS variation on 
timescales less than or of order the Hubble time, traditionally one phrases constraints in terms 
of a constant w. We reproduce this model from UnionOS to serve as a point of comparison. 
Also see UnionOS for models using the standard time varying EOS w{a) = + "^0(1 — a), 
where a = 1/(1 + 2;) is the scale factor, and models with w{z) given in redshift bins. 

In the constant w case the Hubble expansion parameter H = a/a is given by 

H\z)/Hl = fi^(i + zf + + zf^^"-^ + ^^fc(l + zf, (1) 

where VLm is the present matter density, Vt^ the present dark energy density, and Vtk = 
1 — Qm ~ the effective energy density for spatial curvature. 

Figure [T] shows the confidence contours in the w-Qm plane both without and with (min- 
imized in the likelihood fit) spatial curvature. Note that allowing for spatial curvature does 
not strongly degrade the constraints. This is due to the strong complementarity of SN, 
CMB, and BAO data, combined with the restriction to a constant w model. As shown in 
UnionOS, the constraint on curvature in this model is Qk = —0.010 ± 0.012. See UnionOS for 
more plots showing the individual probe constraints. 



4. Braneworld Gravity 



Rather than from a new physical energy density, cosmic acceleration could be due to a 
modification of the Friedmann exp ansion equa.tions arising from an exteii sion of gravitational 
theory. In braneworld cosmology (IDvali et al.ll2000l : iDeffayet et al.ll2002l ). the acceleration is 
caused by a weakening of gravity over distances near the Hubble scale due to leaking into 
an extra dimensional bulk from our four dimensional brane. Thus a physical dark energy is 
replaced by an infrared modification of gravity. For DGP braneworld gravity, the Hubble 
expansion is given by 



Here the present effective braneworld energy density is 



(flat) . 



bw 



, 1 ^m) 



(flat) 



(2) 
(3) 

(4) 
(5) 
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Fig. 1. — 68.3%, 95.4%, and 99.7% confidence level contours on a constant EOS w and 
the matter density fi^ for the individual and combined data sets. The left panel shows 
individual and combined probes in the fiat universe case; the right panel repeats the combined 
systematics contour, and also compares to the statistical only contour, and to the systematics 
contour when simultaneously fitting for spatial curvature. 



and is related to the five dimensional crossover scale = Mp[/(2M|) by fib„ = l/(4ifQr^). 
Note that the only cosmological parameters for this model are VL^ and Vt^ (or fibw), so it has 
the same number of parameters as ACDM. 



The effective dark energy equation of state is given by the simple expression 



w{z) 



1 + ^lm(z) - ^kiz) 



(6) 



where n^iz) = a^{l + zf/{H^/Hl) and nk{z) = Vikil + zY/{H^/Hl). Thus the dark 
energy equation of state at present, wq, is determined by Vtm and VL^] while time varying, it 
is not an independent parameter. So rather than plotting wq vs. Vim or showing constraints 
on the somewhat nonintuitive para r neters Tc or fibw (but see the clear discussion and plots in 
Davis et al.l (120071 ): iRydbeck et al.l (120071 ). though without systematics). Figure [2] illustrates 
the confidence contours in the VLk-VLm plane. This makes it particularly easy to see how 
deviations from flatness pull the value of the matter density. In this and following figures, 
dotted contours show the BAO constraints, dashed for CMB constraints, dot-dashed for SN 
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with systematics, and solid contours give the joint constraints. 



C3 0.0 " 




Fig. 2. — The extradimensional DGP braneworld gravity model does not achieve an accept- 
able fit to the combined data, even allowing for a spatial curvature parameter. The joint 
best fit is in fact a nearly flat model, but with poor goodness of fit: Ax^ = 2.7 relative 
to the ACDM case; also shown is the statistical error only SN contour, which gives a joint 
Ax^ = 15 relative to ACDM. 



For a fiat universe, in order for w to approach —1 the matter density is forced to small 
values. Alternately, pushing the curvature density Qk negative, i.e. introducing a positive 
spatial curvature k, allows w ^ —1 with higher matter density. For a given wq, the amount 
of curvature needed can be derived from Eq. ([6]) to be approximately AQk ~ — Afim/^m, so 
to move a fiat, Q^n = 0.2 universe to Qm = 0.3 requires Qk = —0.5, in agreement with the 
SN contour (being most sensitive to Wq) of Figure [21 

Note that the curvature density cannot exceed 1 — Qm, corresponding to an infinite 
crossover scale Vc, so the likelihood contours are cut off at this line and the region beyond 
is unphysical. However, this does not affect the joint contours. The BAO data contours 
do extend to the limit fi/,. = 1 — Qm] here f2bw = 0, equivalent to the simple OCDM open, 
nonaccelerating universe. 

Most importantly, the three probes do not reach concordance on a given cosmological 
model. The areas of intersection of any pair are distinct from other pairs, indicating that 
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the full data disfavors the braneworld model, even with curvature. This is further quantified 
by the poor goodness of fit to the data, with Ax^ = 2.7 relative to the flat ACDM model 
possessing one fewer parameter, or Ax^ = 4.0 relative to ACDM allowing curvature. This 
indicates the crucial importance of crosschecking probes. Moreover, if we had used only 
the statistical estimates of uncertainties (see the "SN stat" 68% cl contour of Fig. [2j), we 
would have found that Ax^ = 15 rather than 2.7, and possibly drawn exaggerated physical 
conclusions - considering the DGP model 2000 times less likely than it really is, as an 
illustratioE0. Inclusion of systematics is essential for robust interpretation of results. 



5. Doomsday Model 



Perhaps the sim plest g e nerali zation of the cosmological c onstant is the line ar potential 
model, pioneered by iLindd (119861 ) and discussed recently by IWeinberg] (120081 ) . motivated 
from high energy physics. Interestingly, while this gives a current accelerating epoch, in 
the future the potential becomes negative and not only deceleration of the expansion but 
collapse of the universe ensues. Hence the name of doomsday model. 

The potential has two parameters: the amplitude and slope. The amplitude Vq es- 
sentially gives the dark energy density, which is fixed by Qm in a flat universe. (For the 
remainder of the paper we assume a flat universe, for the reasons discussed in ^) The slope 
V = dV/dcj) can be translated into the present equation of state val ue Wq. Thus this is a 



one parameter model in our categorization. See iKal 



cosmological properties of the linear p otential. iLinde ( 19861 ) for a view of it as a perturbation 



osh et al.l (120031 ) for discussion of the 



about zero cosmological constant, and iDimopoulosi (120031 ) for links to the large kinetic term 
ap proach in particl e physics. More recently, this has been considered as a textbook case 
by IWeinbergl (12008! ) . so we will ex amine this model in som e detail. Such dark energy is an 
example of a thawing scalar field (ICaldwell &: Linderll2005l ). starting with w{z ^ 1) = — 1 
and slowly rolling to attain less negative values of w; that is, it departs from A. If it has 
not evolved too far from —1 then its behavior is well described by Wa ~ —1.5(1 + Wq) where 
w{a) = Wo + Wa(l — a). However we solve the scalar field equation of motion exactly for all 
results quoted here. 

As the scalar field rolls to small values of the potential the expansion stops accelerating, 
and when it reaches V = Q then w = 1. However it crosses through zero to negative 



^This is not quite fair as the braneworld model and ACDM model have distinct parameter spaces and 
the reduced x^/ dof is only 1.07 for the statistics only braneworld case. This is one area where Bayesian 
evidence methods, with careful use of priors, would be useful for model comparison. 
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values of the potential, further increasing w, and eventually the dark energy density itself 
crosses through zero, causing w to go to positive and then negative infinity. Thereafter the 
negative dark energy density, acting now as an attractive gravitational force, causes not only 
deceleration but forces the universe to start contracting. The rapid collapse of the universe 
ends in a Big Crunch, or cosmic doomsday in a finite time. 



In the notation used in IWeinbergj (120081 ) . = Vq + (0 — 0o) VJj', with Vq the potential 
energy during the initial frozen state (during high Hubble drag at high redshift) and Vq is 
the constant potential slope. Figure [3] shows the constraints in this high energy physics plane 
Vq-V^'. Note the tight constraints on the initial potential energy Vq, given in units of the 
present critical density. The cosmological constant corresponds to the limit of VJ{ = 0, but 
the slope must always be less than or of order lO"^^'' in Planck units, i.e. unity when shown 
in terms of the present energy density, to match the data. 



2.0 



1.5 



> 
1.0 



0.5 



0.0 
0.0 




Fig. 3. — Constraints on the linear potential model in terms of the high energy physics 
quantities of the primordial amplitude and slope of the potential. Note there is less com- 
plementarity between some of the probes than for the constant w model. Fig. H] translates 
these constraints into ones on the cosmological parameters. 



We can also translate these high energy physics parameters into the recent universe 
quantities of the matter density flm and the present equation of state Wq. Moreover, this is 
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directly related to the doomsday time tdoom, or future time until collapse. A useful approxi- 
mation (though we employ the exact solution) between tdoom, Wq, and the approximate time 
variation Wa = —1.5(1 + Wq) is 

tdoom ~ 0.5H^\1 + Wo)-"-' ^ 0.6H^\-Wa)-''-\ (7) 

Figure H] shows the likelihood contours in the tdoom-^m and Wo-Qm planes. The 95% confi- 
dence limit on tdoom from present observations is 1.24 if q"^, i.e. we are 95% likely to have at 
least 17 billion more years before doomsday! 




0.0 0.1 0.2 ^ 0.3 0.4 0.5 0.0 0.1 0.2 ^ 0.3 0.4 0.5 



Fig. 4. — The future expansion history in the linear potential model has a collapse, or cosmic 
doomsday, at a finite time in the future. The left panel shows the confidence contours for 
the time remaining until collapse; the likelihood contours extend to infinity, with tdoom = oo 
corresponding to the A model. The contours can also be viewed in the equivalent Wo-Qm 
plane (right panel). Current data constraints indicate cosmic doomsday will occur no sooner 
than ~1.24 Hubble times from now at 95% confidence. 



6. Mirage Model 

Given their limited sensitivity to the dynamics of dark energy, current data can appear 
to see a cosmological constant even in the presence of time variation. This is called the 
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"mirage of A" , and we consider mirage models, with a form motivated by the observations 
as discussed below, specifically to test whether the concordance cosmology truly narrows in 
on the cosmological constant as the dark energy. 

Since cosmological distances involve an integral over the energy density of components, 
which in turn are integrals over the equation of state as a function of redshift, there exists 
a chain of dependences between these quantities. Fixing a distance, such as di^s to the 
CMB last scattering surface, can generally lead to an "attractor" behavior in the equation 
of stat e to a common averaged value or the value at a particular redshift. Specifically, 



Linderl (120071 ) pointed out that if CMB data for di^s is well fit by the ACDM model then 
this forces w{z ~ 0.4) ^ —1 for quite general monotonic EOS. So even dark energy models 
with substantial time variation could thus appear to behave like the cosmological constant 
at z ~ 0.4, near the pivot redshift of current data. 

Since current experiments insensitive to time variation inherently interpret the data 
in terms of a constant w given by the EOS value at the pivot redshift, this in turn thus 
leads to the "mirage of A" : thinking that w = — 1 e verywhere, despite models very different 



from A being good fits. See §5.2 of iLinderl (j2008bl ) for further discussion. (Also note that 



attempting to constrain the EOS by combining the CMB diss with a precision determination 
of the Hubble constant Hq only tightens the uncertainty on the pivot equation of state value 
(already taken to be nearly —1) and so similarly does not reveal the true nature of dark 
energy.) 

We test this with a family of "mirage" models motivated by the reduced distance to 
CMB last scattering diss- These correspond to the one parameter subset of the two parameter 
EOS model w{a) = wq + Wa{l — a) with Wa determined by Wa = —3.63(1 + wq). They are 
not exactly equivalent to imposing a CMB prior since diss will still change with Q^n] that is, 
they essentially test the uniqueness of the current concordance model for cosmology: ACDM 
with = 0.28. 

For any model well approximated by a relation Wa = —^(1 + wq), as this model (and 
the previous one) is, the Hubble parameter is given by 

H^H^ = n„ (1 + z)=^ + (1 - r]„) (1 + 2)=^(i+'"«+'"«) e-^"'""/(^+") (8) 

= fi™ (1 + Zf + (1 - Qm) (1 + ^)3(1+-0){1-A) g3A(l+.,o)./(l+.) _ ^g) 



Figure [5] shows constraints in the WQ-Qm plane. It is important to note that w is not 
constant in this model. A significant range of Wq (and hence a larger range of Wa too, roughly 
+0.55 to —1.1 at 68% cl) is allowed by the data, even though these models all look in an 
averaged sense like a cosmological constant. Thus experiments sensitive to the time variation 
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Wa (e.g. a{wa) < 0.36 to know that w{z) is really, not just apparently, within 10% of — 1) 
are required to determine whether the mirage is reality or not. 




Fig. 5. — The mirage subclass of time varying dark energy looks like A in an averaged 
sense. Note that CMB contours are almost vertical, indicating both that the mirage holds, 
preserving the ACDM distance to last scattering, and yet imposes little constraint on wq, 
and hence Wa- Thus the appearance of A does not actually exclude time variation. The 
mirage is broken when the equation of state at high redshift exceeds the matter domination 
value of zero; this causes the wall in the likelihood at wq = A/{1 — A) —1.4. 



7. Vacuum Metamorphosis 



An interesting model where the cosmic acceleration is due to a change in the behavior of 
physical laws, rather than a new physical ene rgy density, is the vacuum metamorph osis model 
([Parker &: Ravall l2000l : ICaldwell et al.ll2006l ). As in Sakharov's induced gravity (ISakharov 
19681 ). quantum fluctuations of a massive scalar field give rise to a phase transition in gravity 
when the Ricci scalar curvature R becomes of order the mass squared of the field, and freezes 
R there. This model is interesting in terms of its physical origin and nearly first principles 
derivation, and further because it is an example of a well behaved phantom field, with 
w < -1. 
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The criticality condition 

R = 6{H + 2H^) = (10) 
after the phase transition at redshift Zf leads to a Hubble parameter 

1 — 

H^/Hl = l]^(l + ^)3 + ___^, z>z,. (12) 

There is one parameter, Vt^, = Qm{zt), in addition to the present matter density flm, 
where 1 — fi^ is proportional to the cosmological constant. The variables Zt and m are 
given in terms of ilm, by Zt = (m^n^/[3f2m(4 — 3^2*)])^/^ — 1 and = 3fim[(4 — 
3r2^)/f2j^/^[(4/m^) — (1/3)]"^/"^. The original version of the model had fixed fi* = 1, i.e. 
no cosmological constant, but if the scalar field has nonzero expectation value (which is not 
required for the induced gravity phase transition) then there will be a cosmological constant, 
and Q-^ deviates from unity. 

Figure [6] shows the confidence contours in the Q^^-Qm plane. To consider constraints 
on the original vacuum metamorphosis model, without an extra cosmological constant, slice 
across the likelihood contours at the fi^ = 1 line. We see that the three probes are inconsis- 
tent with each other in this case, with disjoint contours (indeed the Ax^ = 28.5 relative to 
flat ACDM). Allowing for a cosmological constant, i.e. ^ 1, brings the probes into con- 
cordance, and the best joint flt approaches the lower bound of the region Q^, > Qm- The con- 
dition Q^, = Qrn corresponds to the standard cosmological constant case, with = 1 — ^m, 
since the phase transition then only occurs at Zt = 0. Thus the data do not favor any vacuum 
phase transition. Although this model comprises very different physics, and allows phantom 
behavior, the data still are consistent with the cosmological constant. 



8. Geometric Dark Energy 



According to the Equivalence Principle, acceleration is manifest in the curvature of 
spacetime, so it is interesting to consider geometric dark energy, the idea that the acceleration 
arises from some property of the spacetime geometry. One example of this involves the 
holographic principle of quantum fleld theory as applied to cosmology. This limits the number 
of modes available to the vacuu m energy and so could have an impact on the cosmological 
constant problem (jBoussdl2002l ). The basic idea is that there is a spacelike, two dimensional 
surface on which all the field information is holographically encoded, and the covariant 
entropy bound relates the area of this surface to the maximum mode energy allowed (UV 
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Fig. 6. — The vacuum metamorphosis model involves a phase transition in gravitational laws 
due to quantum effects. Where the quantum field inducing the gravitational deviation has 
no additional zeropoint energy, i.e. cosmological constant, then fi^ = 1, and the data gives 
discordant results. As the model approaches the fi^ = Qm line of pure cosmological constant 
plus matter without a phase transition in the past, the data provide an increasingly good 
fit. (Below the line, the transition takes place further into the future, with no effect on the 
data likelihood.) 



cutoff). The vacuum energy density resulting from summing over modes ends up being 
proportional to the area, or inverse square of the characteristic length scale. However, what 
is perhaps the natural surface to choose, the causal event horizon, does not lead to an energy 
density with accelerating properties. 

Many of the attempts in the literature to overcome this have grown increasingly distant 
from the original concept of holography, though they often retain the name. It is important 
to realize that, dimensionally, any energy density, including the vacuum energy density, has 
p ~ L~^, so merely choosing some length L does not imply any connection to quantum 
holography. We therefore do not consider these models but turn instead to the spacetime 
curvature. 
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8.1. Ricci dark energy R] 



low 



A different approach involves the spacetime curvature directly, as measured through the 
Ricci scalar. This is similar in motivation to the vacuum metamorphosis model of ^ Here 
we consider it purely geometrically, with the key physical quantity being the reduced scalar 
spacetinie cu rvature, in terms of the Ricci scalar and Hubble parameter, as in the model of 
Lindeil booi ). 



n 



R 



ro + ri(l - a). 



(13) 



This quantity directly involves the acceleration. Moreover, we can treat it purely kinemat- 
ically, as in the last equality above, assuming no field equat ions or dynamics. Of course, 
any functional form contains an implicit dynamics (see, e.g., iLinderl (j2008bl )). but we have 
chosen effectively a Taylor expansion in the scale factor a, valid for any dynamics for small 
deviations 1 — a from the present, i.e. the low redshift or low scalar curvature regime. 

At high redshift, as 1 — a is no longer small, we match it onto an asymptotic matter 
dominated behavior for a < = 1 — (1 — 4ro)/(4ri). Solving for the Hubble parameter. 



HVHl 



^4{ro+ri-l) g4ri(l-a) 



a > at 



i I'm Uj , 

The matching condition determines 

'4ro + 4ri — 



a < at . 



4ro+4ri — 1 



4ri 



,l-4ro 



(14) 
(15) 



(16) 



so there is only one parameter independent of the matter density. 

Note also that we can define an effective dark energy as that part of the Hubble pa- 
rameter deviating from the usual matter behavior, with equation of state generally given 

by 



47^ 



wia] 



- J;^ida/a){l-4n) 



(17) 



For the particular form of Eq. f[T3|) we have 



Wo 



1 -4ro 

3(1 



This model has one EOS parameter in addition to the matter density. We can therefore 
explore constraints either in the general kinematic plane vq-ti, or view them in the Qm-Wo 
plane. Figure [7| shows both. 
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Fig. 7. — Geometric dark energy in the Riow niodel describes the acceleration directly through 
the reduced Ricci scalar, or spacetime curvature. This can be viewed in a kinematic sense, in 
the tq-ti plane, or in a dark energy sense in the Qm-Wo plane. The data favor wq = —1 but 
this is not A, instead representing distinct physics. For rg+ri > 1/4, above the diagonal line, 
early matter domination is violated, and the CMB and BAO likelihoods avoid this region, 
as seen in the left panel; the matter density also cannot then be uniquely defined so the 
equivalent region is excluded from the right panel. 

Good complementarity, as well as concordance, exists among the probes in the tq-ti 
plane. One obtains an excellent fit with (ro,ri) = (0.81, —0.72). The value of 71 today, tq, 
approaches unity, the deSitter value. Recall that 71 = 1/4 corresponds to matter domina- 
tion, and 7?. = 1/2 to the division between decelerating and accelerating expansion, so this 
kinematic approach clearly indicates the current acceleration of the universe. 

An interesting point to note is that ACDM is not a subset of this ansatz, i.e. the physics 
is distinct. No values of tq and ri give a ACDM cosmology. However, the Hubble diagram 
for the best fit agrees with that for ACDM to within 0.006 mag out to z = 2 and 0.3% in 
the reduced distance to CMB last scattering. This is especially interesting as this geometric 
dark energy model is almost purely kinematic. The agreement appears in the Qm-wo plane as 
contours tightly concentrated around Wq = —1, despite there being no actual scalar field or 
cosmological constant. Again we note the excellent complementarity between the individual 
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probes, even in this very different model. 



8.2. Ricci dark energy Rhigh 

Rather than expanding the spacetime curvature around the present value we can also 
consider the deviation from a high redshift matter dominated era. That is, we start with 
a standard early universe and ask how the data favors acceleration coming about. In this 
second geometric dark energy model (call it Rhigh for high redshift or large values of scalar 
curvature), the value of TZ evolves from 1/4 at high redshift. From the definition of TZ, it 
must behave asymptotically as 

5H^' 



4 



1 — 3Wf, 



i [1 + 4aa-2'"-], (19) 



where 5H^ = {H"^ / Hi)— Vim {1+zY is the deviation from matter dominated behavior, and Woo 
is the associated, effective equation of state at high redshift, approximated as asymptotically 
constant. 

Next we extend this behavior to a form that takes the reduced scalar curvature to a 
constant in the far future (as it must if the EOS of the dominant component goes to an 
asymptotic value): 

_ 1 aa~^'"°° , , 

So today IZ = l/4 + a/(l + /?) and in the future IZ = 1/4, + a/ [3. By requiring the correct 
form for the high redshift Hubble expansion, one can relate the parameters a and j3 by 

a = (3/5ti;oo/4) [In ln(l + /5)] , (21) 

and finally 

H^/Hl = n^a-^ (1 + /3a-3"'-)-i^f^-/i'^(i+^). (22) 

The Rhigh geometric dark energy model has two parameters /? and Woo, in addition to the 
matter density flrn- This is the first such model we consider, and all remaining models also 
have two EOS parameters. Although current data cannot in general satisfactorily constrain 
two parameters, and so for all remaining models we do not show individual probe constraints, 
if the EOS phase space behavior of the model is sufficiently restrictive then reasonable joint 
constraints may result. 

Figure [8] shows the joint likelihoods in the Qm-Woo and flm-P planes, with the third 
parameter minimized over. We see that the data are consistent with the cosmological con- 
stant behavior w^o = — 1 in the past (this is only a necessary, not sufficient condition for 
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ACDM), and indeed constrain the asymptotic high redshift behavior reasonably well, in par- 
ticular to negative values of Woo- This indicates that the Ricci scalar curvature definitely 
prefers a nearly-standard early matter dominated era, i.e. the deviations faded away into the 
past. This has important implications as well for scalar-tensor theories that would modify 
the early expansion hist ory; in particula r , the data indicate that deviations in TZ must go 



approximately as (see iLinder fc CahnI (120071 )) not sometimes assumed. 



The parameter f3 helps determine the rapidity of the Ricci scalar transition away from 
matter domination. This varies between f3 = 0, giving a slow transition but one reaching 
a deceleration parameter q = —oo in the asymptotic future, and 3> 1, giving a rapid 
deviation but with smaller magnitude. A cosmological constant behavior has /? ~ 3, as 
discussed below. 




Fig. 8. — Geometric dark energy in the Rhigh model describes the acceleration directly 
through the reduced Ricci scalar curvature and deviations from early matter domination. 
The left panel shows the Qm-Woo plane, indicating the nature of the deviation (wqo = 
corresponds to no transition away from matter domination), and the right panel shows the 
^m-/3 plane, indicating the rapidity and fate of the deviation. The curve in the left panel 
corresponds to whether the fate of the universe is de Sitter; we also show the individual 
probe constraints, fixing (3 to the de Sitter value (not minimizing over j3 as for the joint 
contour), to show that SN closely map the fate of the universe. In the right panel the curve 
is the cut through parameter space, fixing Woo = — 1, corresponding to ACDM. 



- 20 - 



Within the three dimensional parameter space, two subspaces are of special interest. 
One is where Woo = —1, a necessary condition for consistency with A, as mentioned. The 
other corresponds to a deSitter asymptotic future, defined by the line 

/?dcs = - 1. (23) 

Note that unlike the previous geometric model Riow, the Rhigh model does include A as the 
limit when both these conditions are satisfied, Woo = — 1 and (3 = — 1. This relation for 
the A limit is shown as a curve in the 1^^-/3 plane. There is an overlap with the joint data 
likelihood, though one must be careful since the contours have been minimized over Woo- 

Interestingly, we can actually use the data to test consistency with a de Sitter asymptotic 
future. This is shown by the curve in the flm-Woo plane. We see that SN are the probe most 
sensitive to testing the fate of the universe, with the SN contour oriented similarly to the 
curve given by Eq. ( l23i) that passes through the best fit. Thus the data are consistent 
with Woo = — 1 and with a de Sitter fate separately, though some tension exists between 
satisfying them simultaneously. Thus, this geometric dark energy may be distinct from the 
cosmological constant. 



9. PNGB Model 



Returning to high energy physics models for dark energy, one of the key puzzles is how 
to prevent quantum corrections from adding a Planck energy scale cosmological constant or 
affecting the shape of the potential. This is referred to as the issue of technical naturalness. 
Pseudo-Nambu Goldstone boson (PNGB) models are technically natural, due to a shift 
symmetry, a nd so can be co nsidered strongly physically motivated (perhaps even more so 
than A). See iFriemarj ( 19951) for an early co s mological analysis of PN GB as dark energy and 
more recent work by iDutta fc Sorbd (120071 ) ; lAbrahamse et al.l (120081 ) . 



The potential for the PNGB model is 



l^(0) = K[l + cos(0//)], 



(24) 



with K setting the magnitude, / the symmetry energy scale or steepness of the potential, 
and (pi is the initial value of the field when it thaws from the high redshift, high Hubble drag, 
frozen state. These three parameters determine, and can be thought of as roughly analogous 
to, the dark energy density, the time variation of the equation of state, and the value of the 
equation of state. The dynamics of this class of models is sometimes approximated by the 
simple form 



wia) 



-1 + (1 + Wo)a^ 



(25) 
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with F roughly inversely related to the symmetry energy scale /, but we employ the exact 
numerical solutions of the field evolution equation. 

PNGB models are an example of thawing dark energy, where the field departs recently 
from its high redshift cosmological constant behavior, evolving toward a less negative equa- 
tion of state. Since the EOS only deviates recently from w = —1, the precision in measuring 
wq is more important than the precision in measuring an averaged or pivot EOS value. 
SN data provide the tightest constraint on wq. In the future the field oscillates around its 
minimum with zero potential and ceases to accelerate the expansion, acting instead like 
nonrelativistic matter. 

Figure [9] illustrates the constraints in both the particle physics and cosmological param- 
eters. The symmetry energy scale could provide a key clue for revealing the fundamental 
physics behind dark energy, and it is interesting to note that these astrophysical observa- 
tions essentially probe the Planck scale. For values of / below unity (the reduced Planck 
scale), the potential is steeper, causing greater evolution away from the cosmological con- 
stant state. However, the field may be frozen until recently and then quickly proceed down 
the steep slope, allowing values of Wq far from —1 but looking in an average or constant w 
sense like (w) ~ —1. Small values of (pi/f have the field set initially ever more finely near 
the top of the potential; starting from such a fiat region the field rolls very little and w stays 
near —1 even today. In the limit (pi/f = the field stays at the maximum, looking exactly 
like a cosmological constant. The two effects of the steepness and initial position mean that 
the cosmological parameter likelihood can accommodate both wo ~ — 1 and wq approaching 
as consistent with current data. However, to agree with data and 1 + ~ 1 requires 
/ <^ 1 and fine tuning - e.g. for / = 0.1 one must balance the field to within one part in a 
thousand of the top. Thus in the left panel there exists an invisibly narrow tail extending 
along the y-axis to / = 0. In the right panel, we show how taking more natural values 
/ > 0.5 removes the more extreme values of Wq caused by the unnatural fine tuning. 



10. Algebraic Thawing Model 



While PNGB models involve a pseudoscalar thawing field, we can also consider scalar 
fields with thawing behavior. Any such fields that are neither fine tuned nor have overly steep 
potentials must initially depart from the cosmological constant behavior along a specific track 
in the EOS ph ase space, characterized by a forn i of slow roll behavior in the matter domi- 



nated era. fSee lCaldwell fc Linded fl2005f ): lLindei] (120061 ) :IScherrer fc Sen I (120081 ) : ICahn et al. 



(120081 ).) Here we adopt the algebraic thawing model of iLinderl (l2008al ). specifically designed 
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Fig. 9. — Left panel: PNGB model dynamics involves a competition between the steepness 
of the potential, given by the symmetry energy scale /, and the initial field position If 
the potential is very steep, / -C 1, the field will roll so rapidly to the potential minimum that 
the dark energy density never becomes significant, unless (pi/ f is fine tuned very near zero. 
For natural energy scale values near the Planck scale, / ~ 1, a wide variety of (pi/ f are viable. 
Right panel: The field spends a long period frozen, acting as a cosmological constant before 
thawing and evolving to a present EOS wq. For steep potentials with / <^ 1, the thawing 
can be rapid and result in evolution to Wq far from —1, yet still be consistent with data. 
The solid confidence level contours in the WQ-Vtm plane show PNGB results for energy scales 
/ > 0.1, while the white outline contours consider only PNGB models with more natural 
energy scales / > 0.5; the latter favors models closer to the cosmological constant behavior. 

to incorporate this physical behavior: 



l + w 




(26) 



H^/Hl 



VLm a ^ + (1 — rim) exp 



{1 - (l-a + aa^)^/^} , (27) 



where a = 1/(1 + 6) and b = 0.3 is a fixed constant not a parameter. The two parameters 
are Wq and p and this forr n fulfills the physi cal dynamics condition not only to leading but 
also next-to-leading order ( Cahn et al.lboosi ). 
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The physical behavior of a minimally coupled scalar field evolving from a matter dom- 
inated era would tend to have p G [0,3]. Since we want to test whether the data points 
to such a thawing model, we consider values of p outside this range. Results are shown in 
Figure [101 

For p < 0, the field has already evolved to its least negative value of w and returned 
toward the cosmological constant. The more negative p is, the less negative (closer to 0) the 
extreme value of w is, so these models can be more tightly constrained as p gets more strongly 
negative. As p gets more positive, the field takes longer to thaw, increasing its similarity to 
the cosmological constant until recently, when it rapidly evolves to wq. Such models will be 
very difficult to distinguish from A. If we restrict consideration to the physically expected 
range p G [0,3], this implies Wq < —0.57 at 95% confidence in these thawing models, so 
considerable dynamics remains allowed under current data. This estimation is consistent 
with the two specific thawing models already treated, the doomsday and PNGB cases. 

The goodness of fit to the data is the best of all models considered here, even taking into 
account the number of fit parameters. This may indicate that we should be sure to include 
a cosmological probe sensitive to Wq (not necessarily the pivot EOS Wp) and to recent time 
variation Wa-, such as SN, in our quest to understand the nature of dark energy. 



11. Early Dark Energy 



The other major class of dark energy behavior is that of freezing models, which start out 
dynamical and approach the cosmological constant in their evolution. The tracking subclass 
is interesting from the point of view again of fundamental physics motivation: they can 
ameliorate the fine tuning problem for the amplitude of the dark energy density by having 
an attracto r behavior iii their dynamics, drawing from a large basin of attraction in initial 
conditions ( Zlatev et al.l 1999 ). Such models generically can have nontrivial amounts of dark 
energy at high redshift; particularly interesting are scaling models, or tracers, where the dark 
energy has a fixed fraction of the energy density of the dominant compo nent. These cari be 
motivated by dilatation symmetry in particle physics and string theory (jWetterichlll988l ). 



As a specific model of such early dark energy we adopt that of iDoran fc Robbers I (120061 ) 



with 



DE(a) 



+ ^^e(l 



-3^«o^ 



(28) 

for the dark energy density as a function of scale factor a = + z). Here = 1 — 
is the present dark energy density, fig is the asymptotic early dark energy density, and Wq is 
the present dark energy EOS. In addition to the matter density the two parameters are fie 
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Fig. 10. — Algebraic thawing model incorporates the expected physical behavior of a thawing 
scalar field rolling slowly from a matter dominated era. Such a model is a fairly generic 
parametrization for this class of physics when p G [0, 3], and has a strong goodness of fit. 



and Wo. 

The Hubble parameter is given by H^/Hq = /[^ — ^be{0')]- The standard formula 

for the EOS, w = — 1/(3[1 — i7DE(«)]) d\nQD-E{a)/d In a, does not particularly simplify in this 
model. Note that the dark energy density does not act to accelerate expansion at early times, 
and in fact w —>■ 0. However, although the energy density scales like matter at high redshift, 
it does not appreciably clump and so slows growth of matter density perturbations. We will 
see this effect is crucial in constraining early dark energy. 

Figure [TT] shows the constraints in the Qm-^e and ^Ic-Wq planes. Considerable early dark 
energy density appears to be allowed, but this is only because we used purely geometric in- 
formation, i.e. distances and the acoustic peak scale. The high redshift Hubble parameter for 
a scaling so lution is multiplied b y a factor — He relative to the case without early dark 
energy (see 



Doran et al. This means that the sound horizon is shifted according to 



s ~ a/I — f2e, but a geometric degeneracy exists whereby the acous tic peak angular sc a le can 
be preserved by changing the value of the matter density (see iLinder fc Robbers J2008h 
for a detailed treatment). This degeneracy is clear in the left panel. 
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However, as mentioned, the growth of perturbations is strongly affected by the unclus- 
tered early dark energy. This suppresses growth at early times, leading to a lower mass 
amplitude cxg today. To explore the influence of growth constraints, we investigate adding a 
growth prior of 10% to the data, i.e. we require the total linear growth (or ag) to lie within 
10% of the concordance model. The innermost, white contour of the left panel of Fig. [TT] 
shows the constraint with the growth prior. In the right panel we zoom in, and show Qe 
vs. Wq, seeing that the degeneracy is effectively broken. The amount of early dark energy is 
limited to fi p < 0.03 8 at 95% cl. Similar conclusions were found in a detailed treatment by 
Doran et al.l J2007bh . 



We find a convenient fitting formula is that for an early dark energy model the total 
linear growth to the present is suppressed by 



9o 



0.01 



X 5.1%, 



(29) 



relative to a model with fig = but all other parameters fixed. Thus appreciable amounts 
of early dark energy have significant effects on mat ter perturbations, a nd w e might expect 



nonlinear growth to be even more sensitive (e.g. see iBartelmann et al.l (120061 ) ). 



12. Growing Neutrino Model 

While freezing or scaling models such as the early dark energy model just considered are 
interesting from the physics perspective, they generically have difficulty in e volving naturally 



to sufficiently neg a tive E OS by the present. The growing neutrino model of lAmendola et al. 



( 120071 ) : IWetterichI (120071 ) solves this by coupling the scalar field to massive neutrinos, forcing 
the scalar field to a near cosmological constant behavior when the neutrinos go nonrela- 
tivistic. This is an intriguing model that solves the coincidence problem through cosmo- 
logical selection (the time when neutrinos become nonrelativistic) rather than tuning the 
Lagrangian. 

The combined dark sector (cosmon scalar field plus mass-running neutrinos) energy 
density is 

o / A fidsa3 + 2fi,(a3/2-a3) 

^^^^"^ = l-fids(l-a3) + 2fi.(a3/2-a3)' « > (^0) 
fids (a) = fie, a<at, (31) 

where fids = 1 — fim is the present dark sector energy density. The Hubble parameter can 
be found by H^/Hq = fim,o^^/[l — fids(fl)] as usual. The two free dark parameters are the 
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Fig. 11. — Early dark energy represents an example of a freezing model with interesting 
particle physics motivations. The left panel shows the constraints on Qe and flm from purely 
geometric data, as used throughout this article. The degeneracy evident in the contours 
leaves the acoustic scales unchanged, but hides the shift in the sound horizon caused by 
early dark energy, leading to possible misinterpretation of the correct cosmological model. 
The degeneracy can be broken by adding growth information, here a 10% prior on total 
linear growth (or ag), as shown by the white outline contours. This tightly restricts the 
early dark energy density to contribute no more than a few percent. The right panel shows 
the ^le-Wo constraints including the growth prior. 



neutrino mass or density Q,y = m,y{z = 0)/(30.8/i^ eV) and the early dark energy density Qe- 
The transition scale factor at is determined by intersection of the two behaviors given for 
fids (a)- 

The equation of state is 

with w = before the transition, i.e. a return to the standard early dark energy model. One 
can therefore translate fl^ or m^{z = 0) into Wq = —1 + fi^^/fids = — 1 + ^u/ (1 — ^m)- 

Figure [12] shows the constraints in the m,y{z = 0)-Qe plane. As in the previous early 
dark energy model, the geometric degeneracy is clear. Again, when we add growth infor- 
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mation in the form of a 10% prior on the total hnear growth (or the mass variance erg), 
the constraints tighten considerably, as shown in the right panel. The 95% confidence level 
limit on the neutrino mass from this current cosmological data is then 2.1 {h/O.?)"^ eV (1.2 if 
only statistical uncertainties are taken into account). These limits are comparable to astro- 
physical constraints from similar types of data applied to standard, constant mass neutrinos 
( iGoobar et al.ll2006l : iTegmark et al.ll2006l ). Note that because the neutrino mass grows due 
to the coupling, the value toda y can actually be la rger than that at, say, z ~ 3 where Lyman 
alpha forest constraints apply (ISeljak et al.l 120061 ). 




[(h/0.7)"eV] 



Fig. 12. — Growing neutrino model, coupling a dark energy scalar field to massive neutrinos, 
can solve the coincidence problem. The left panel shows the constraints from purely geomet- 
ric data, while the right panel (note the different vertical scale) adds a 10% prior on total 
linear growth (also see Fig. [TTl) . The neutrino mass today becomes tightly constrained to an 
interesting range, and comparison with laboratory limits could lead to evidence of varying 
neutrino mass. 



13. Conclusion 

We have considered a wide variety of dark energy physics quite different from the cos- 
mological constant. These include a diversity of physical origins for the acceleration of the 
expansion: from dynamical scalar fields to dark energy that will eventually cause deceleration 
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and collapse, to gravitational modifications arising from extra dimensions or from quantum 
phase transitions, to geometric or kinematic parametrization of the acceleration, to dark 
energy that may have infiuenced the early universe and that may have its magnitude set by 
the neutrino mass. The comparison to ACDM and constant w cases covers 5 one-par ameter 



and 5 two-parameter dark energy equation of state models. (ILinder fc Hutererl (120051 ) detail 
how even next generation data will not generically be able to tightly constrain more than 
two such parameters.) 

Two key results to emphasize are that current data 1) are consistent with A, and 2) are 
also consistent with a diversity of other models and theories, even when we restrict consider- 
ation to those with at least modest physical motivation or justification. As explicitly shown 
by the mirage model, any inclination toward declaring A the answer based on consideration 
of a constant w has an overly restricted view. The need for next generation observations 
with far greater accuracy, and the development of precision growth probes, such as weak 
gravitational lensing, is clear. All major classes of physics to explain the nature of dark 
energy are still in play. 

However there are already quite hopeful signs of imminent progress in understanding the 
nature of dark energy. For example, for the braneworld model tight control of systematics 
would decrease the goodness of fit to Ax^ = +15, even allowing for spatial curvature, 
diminishing its likelihood by a factor 2000 naively, effectively ruling out the model. For 
the doomsday model, improving errors by 30% extends our "safety margin" against cosmic 
collapse by 10 billion years - a nonnegligible amount! Every improvement in uncertainties 
pushes the limits on the neutrino mass within the growing neutrino model closer toward other 
astrophysical constraints - plus this model essentially guarantees a deviation from w = —1 
of 0.1 {m^Q/eV), excitingly tractable. Terrestrial neutrino oscillation bounds already provide 
within this model that 1 +w > 0.005. 

As points of interest, we note that the model with noticeably positive Ax^ relative 
to A, and hence disfavored, is completely distinct from the cosmological constant, i.e. the 
braneworld model has no limit within its parameter spaces equivalent to A. This does not 
say that no such model could fit the data - the i?iow niodel is also distinct from A but 
fits as well as many models. Certainly many successful models under current data do look 
in some averaged sense like a vacuum energy but this does not necessarily point to static 
dark energy. Two serious motivations to continue looking for deviations are that physicists 
have failed for 90 years to explain the magnitude required for a cosmological constant, and 
that the previous known occurrence of cosmic acceleration - infiation - evidently involved a 
dynamical field not a cosmological constant. 

To guide further exploration of the possible physics, we highlight those models which 
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do better than A: the geometric dark energy and algebraic thawing approaches. One of the 
sole models where adding a degree of freedom is justified (albeit modestly) by the resulting 
reduction in is the Rhigh model directly studying deviations of the spacetime curvature 
from the matter dominated behavior. This has one more parameter than the constant w 
EOS approach, but improves in by 1. In addition, it has a built-in test for the asymptotic 
de Sitter fate of the future expansion. We recommend that this model be considered a 
model of interest for future fits. The other model improving by at least one unit of 
is the algebraic thawing model, performing better than the other thawing models, with a 
general parametrization explicitly incorporating the physical conditions imposed by matter 
domination on the scalar field dynamics. 

The diversity of models also illustrates some properties of the cosmological probes be- 
yond the familiar territory of vanilla ACDM. For example, for the algebraic thawing and 
other such evolutionary models, the premium is on precision of Wq and Wa much more than 
the averaged or pivot EOS value Wp. Not all models possess the wonderful three-fold com- 
plementarity of the probes seen in the constant w case; for many of the examples BAO and 
CMB carry much the same information as each other. However, we clearly see that for every 
model SN play a valuable role, complementary to CMB/BAO, and often carries the most 
important physical information: such as on the doomsday time or the de Sitter fate of the 
universe or the Planck scale nature of the PNGB symmetry breaking. 

The diversity of physical motivations and interpretations of acceptable models highlights 
the issue of assumptions, or priors, on how the dark energy should behave. For example, in 
the i?iow model should priors be flat in tq, ri or in Qm, wo] in the PNGB model should they 
be flat in /, or in Qm, ujq, etc.? Lacking clear physical understanding of the appropriate 
priors restricts the physical meaning of any Bayesian evidence one might calculate to employ 
model selection; the goodness of fit used here does not run into these complications that 
can obscure physical interpretation. 

We can use our diversity of models for an important consistency test of our understand- 
ing of the data. If there would be systematic trends in the data which do not directly project 
into the ACDM parameter space (i.e. look like a shift in those parameters), then one might 
expect that one of the dozen models considered might exhibit a significantly better fit. The 
fact that we do not observe this can be viewed as evidence that the data considered here 
is not flawed by significant hidden systematic uncertainties. The data utilize the UnionOS 
compilation of uniformly analyzed and crosscalibrated Type la supernovae data, constitut- 
ing the world's published set, with systematics treated and characterized through blinded 
controls. The data are publicly available at http://supernova.lbl.gov/Union, and will be 
supplemented as further SN data sets become published; the site contains high resolution 
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figures for tliis paper as well. 

However, to distinguisli deeply among tlie possible physics behind dark energy requires 
major advances in several cosmological probes, enabling strong sensitivity to the time vari- 
ation of the equation of state. This is especially true for those models that are now or were 
in the past close to the cosmological constant behavior. We are getting our first glimpses 
looking beyond A, but await keen improvements in vision before we can say we understand 
the new physics governing our universe. 
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